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GEOMETRIC ERGODICITY OF A HYPOELLIPTIC DIFFUSION 
MODELLING THE MELT-SPINNING PROCESS OF NONWOVEN 

MATERIALS 

MARTIN KOLB, MLADEN SAVOV, AND ACHIM WUBKER 

Abstract. We analyze the large time behavior of a stochastic model for the lay-down of 
fibers on a conveyor belt in the production process of nonwovens. It is shown, that under 
weak conditions this degenerate diffusion process is strong mixing, confirming a conjecture of 
Grothaus and Klar. Moreover, under some additional assumptions even geometric ergodicity 
is established using probabilistic tools - described in the book of Meyn and Tweedie - in 
combination with methods from stochastic analysis. 

o 

P 

OO ! 1. Introduction and Notation 

<N 

Motivated by industrial production methods of nonwoven materials, several mathematical 
models of different complexity have been developed in [1], [11] and [4] in order to accurately 
describe and optimize these web forming processes. In the melt-spinning process a large 
amount of polymer fibers are exposed to an turbulent air flow, which induces an entanglement 
of the different fibers and eventually the formation of a web. Production methods of this type 
are used in the textile, hygiene as well as the building industry. Typical products include 
clothing textiles, industrial filters and insulating materials. For more details concerning the 
engineering background as well as the modeling framework we refer to [10] and references 
therein. In this work we are investigating the simplest of these models. This stochastic 
model for the fiber lay down process has been presented in [4] . It is described by a stochastic 
differential system modeling the image of the fiber motion on the conveyer belt under the 
influence of the turbulent air flow. The interaction of the fiber motion and the air flow is 
^>>0 | formulated via the following stochastic differential equation 

dCt = T(a t ) dt 

da t = adB t -V<p{t t )-T{at) ± dt, 
where r(a) = (cos(a), sin(a)) T and r(a)' L = (— sin(a), cos(a)) . Let us emphasize that it is 
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natural to consider the driving Brownian motion (Bt)t>o an d thus also the process {at)t>o as 
processes on the circle §i or alternatively in [0, 2ir] with periodic boundary conditions. The 
natural state space S for the SDE (1) is therefore S = M. 2 x Si. The generator L v of this 
process (X t ) t >o = (6,at)t>o is of course given by 



X 

5-S 



^2 

L v = -^d 2 a + cos(a)% + sin(a)% - V </?(£) • r(a)~ L <9 Q . 

This differential operator is obviously neither sectorial nor elliptic and therefore the usual 
properties are much more difficult to establish than in the standard elliptic situation. Usually 
even the proof of existence of an invariant distribution and the derivation of its basic properties 
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7 f t f(X s )ds-E,[f] < 1 d + ^HZ-E^/]! 

1 JO L 2 (P a ) V c \ l VI/ 



can be very demanding. In the case of (1) fortunately a straightforward calculation - using 
the generator - (see Proposition 2.4 below as well as Theorem 3.2 in [6]) shows that the 
measure fi on S, given by 

(1.2) (jt(d€, da) := ^e"^ d£ da 

is a promising candidate for an invariant distribution for the process (£ti a t)t>o, once 

N := f e-v^dtda < oo. 
Js 

In the recent paper [6] M. Grothaus and A. Klar use analytic methods close to those presented 
in [19] in order to show that 

IV) 

Actually, Grothaus and Klar are even able to provide a rather explicit description of the 
constants involved. On the one hand this result is obviously quite satisfactory as it gives 
explicit upper bounds on the rate of convergence. On the other hand it only makes an 
assertion about the large time behavior of time averages and not of the process itself. The 
authors of [6] conjecture that for smooth potentials the process (X t )t>o is in fact also strong 
mixing in the sense of 

lim / \E x [f(X t )} - E M [/]| 2 ^(x) = 0. 

We will prove this for a larger class of potentials tp. A rather immediate question is, whether 
one has a geometric rate of convergence. This will be answered affirmatively in section 3 
under further assumptions on tp in the sense that under appropriate conditions Theorem 3.2 
implies 

- lim \ \og\\¥ x {X t G •) -n\\ TV > 

for every x £ S. In contrast to [6] we use mainly probabilistic techniques and do not rely on 
Dirichlet form and semigroup methods. A geometric rate of convergence is usually associated 
to the spectral gap property in some suitable chosen norm (see e.g. [14], [16], [18] and [26]) and 
therefore one might try to prove this using methods from spectral theory. The non-ellipticity 
of this process makes such an approach rather difficult and we do not know, whether recently 
developed analytic tools (see [24]) might be applied in the situation at hand. In any case, 
our approach only relies on rather basic arguments much more elementary than those needed 
in the theory of hypocoercivity. Our approach has the potential to extend to the situation 
of a moving conveyor belt (see e.g. equation (6) in [10]), where the existence of a stationary 
distribution is not known apriori. Last but not least we believe that our stochastic analysis 
shows once more the power of the coupling approach outlined in [16] and offers additional 
insight into the probabilistic mechanism underlying the fiber lay down process. 

The plan for this paper is the following : In the section 2.1 we present existence, uniqueness 
as well as some basic properties of the solution of the solutions to . Some of these are known 
but we give a complete account in order to be independent from [6], where it is made use of 
the theory of generalized Dirichlet forms; others contain new results such as e.g. the strong 
Feller property proved in Proposition 2.3. In section 2.2 we prove ergodicity under weak 
conditions on <p and in section 3 we show geometric ergodicity under stronger assumption on 
(p. 
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2. First Properties of The Fiber Lay Down Process 

We start by establishing existence and uniqueness of solutions of (1) in the weak sense. 

Proposition 2.1. Assume that <p G C 1 (M?). For every initial state x = (£ci; a o) £ S the 
stochastic differential equation (1) has a unique weak solution ¥ x , which is non-explosive. Let 
us denote by the solution for the case <p = then for every t > we have ¥ x \ Tt « 
Q x \ J~t and the following Girsanov formula holds: 

(2-1) ^ \T t = exp(- J V^ s )-r(a s )da s - Y J |V<^(£ S ) ■ r(a s )\ 2 ds) 

Proof. The proof essentially follows ideas outlined in [25] and is devided into two simple steps: 
Step 1: We first look at the situation (p = 0. In this case the stochastic differential equation 
is well known to have (explicit) solutions, which are unique in the weak as well as in the 
pathwise sense. The solution forms a conservative strong Markov process (£,ti a t)t>0 having 
the Feller property in the sense that for every bounded measurable g : S — > M and t > the 
function 



S3x^ J </((#, a?)), 



is continuous. 

Step 2: We now look at the case ip E C 1 (1R 2 ) assuming additionally that V92 is bounded. In 
this situation we rely on the Cameron-Martin-Girsanov formula as stated e.g. as Theorem 9.1 
in [20]. According to this theorem or respectively Corollary 9.3 in [20] we conclude that the 
martingale problem associated to the operator has a unique solution (P x )xes satisfying 

(2.2) F x = expf- J*V<p(&) ■ r(a° s ) da° s - ^ j)^^s) ■ r(a° s )\ 2 daW 

where we denote by Q x the law of the Markov process a® )t>o from Step 1 with the initial 
condition (£o> a o) = x - ^ n particular the exponential in (2.2) is a true martingale. 
Step 3: We now drop the condition of boundedness of Vy. 

Let us denote by 



<tr = inf{t >0\X t = a t ) i B(0, R) x Si} 



the first exit time from the compact set B(0, R) x §1. We now show that for every x £ S and 
t > 

(2.3) lim FJo-r >t)=l. 

Assume for a moment that (2.3) holds true, then 



(2.4) / M t dQ x > / M tAaR dQ x = ¥ x (a R > t) 1 

as R — > 00. This shows that the local martingale (Mt)t>o is a true martingale. 

In order to prove (2.3) we rely on the Liapunov method as presented e.g. in section 6.7 
[20]. Let us set tf(f,a) = |£| 2 + 1, then for all £ = (£1,6) with > 1 

L^(^a) = cos(a)%*(£, a) + sin(a)% tt(£,a) < 2(£ x +£ 2 ) < 2tt(f,a). 
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Consequently, 



defines a P x -supermartingale and therefore 



(2.5) E x [e- 2 ^^(X tAaR )] <*(x). 

Using the fact that ^(X aR ) > R 2 + 1 on {tr < oo} we get 



R 2 + l 



which implies (2.3). 



□ 



Remark 2.2. The fact that the process (Xt)t>o does not explode is more or less obvious due 
to the fact that the 'velocity' in the ^-component is bounded by one. Thus the process can 
not reach infinity in finite time. This is reflected in the elementary choice of the Liapunov 
function and the fact that the a-component does not play any role here. We included the proof 
here as the Liapunov function method will most probably be needed in the analysis of more 
complex models of the fiber lay down process. 

2.1. Feller and Recurrence Properties. We now establish smoothing properties of the 
stochastic process (£t,&t)t>o- These are essential in order to prove the mixing properties in 
the next section. Again our argument is perturbation theoretic in the sense, that we look for 
conditions, which allow to transfer these properties known to hold in the case (p = to the 
general situation. Our arguments are again inspired by the corresponding ones in [25] (see 
also [12] for the use of Girsanov's formula in the derivation of Feller properties). 

Proposition 2.3. Assume that ip G C 1 (IR 2 ) and let (Pt(x,dx'))t>o denote the transition 
probability kernels of the Markov process ((Xt)t>o, (P^zes)- Then for every t > and x £ S 
one has 



is continuous. In particular, Pt is strong Feller for each t > 0. 

Proof. The proof very closely follows Proposition 1.2 in [25]. Consider the case ip = first. 
In this case the generator is obviously hypoelliptic and by Hormander's theorem (see e.g. 
chapter 7 in [23]) we conclude that there exists q% € C°°(S x S) such that the transition 
probability Qt(x,dx') of Q satisfies 



Since x — > Q x is continuous with respect to the weak convergence of probability measures we 
conclude the strong Feller property of Qt for every t > 0. 

Now let (p satisfy the conditions of Proposition (2.1). Then according to the Girsanov 
formula established in the proof of this result we conclude that there exists pt(x, x') such that 
Pt(x, dx') = pt(x, x')dx' for every fixed t > 0. As in [25] we conclude using a result of Revuz 
that it suffices to show, that P t f £ Cb(S) for every bounded and measurable / : S — > R and 



Pt(x,dx') = pt(x, x')dx' , pt(x,x) > 



and 



SBx^p t {x,-) £L'(S) 



Qt(x,dx') = qt(x,x')dx' 



t > 0. 
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Let (Bt)t>o be a Brownian motion in [0, 2tt] with periodic boundary conditions defined on 
the filtered probability space (fi, J 7 , (J r t)t>QiQ)- We denote by (X® (x)) s >o the strong solution 
of (1) with ip = and with initial condition Xq = x, i.e. 

6 = . ' ' and a t = B t . 
VJo sin(-B s ) dsj 

Then we have 

(2.6) P t f(x)=E ¥ °[f(X t )] = E®[f(X?(x))M t (x)], ^ 

where Mt(x) is the exponential martingale given by the Girsanov formula in Proposition 2.1, 
but with ((,s,O(s)s>0 substituted by the (Xg(x)) s >o- Now let (x n ) nS N C S denote a sequence 
converging to x £ S. Then we have (Xg(x n )) s >o — > (X®(x)) s >o locally uniformly in Im- 
probability and by convergence properties of the stochastic integral we have Mt{x n ) — > Mt(z) 
in Q-probability. As for n G N the processes {Mt{z n ))t and {Mt{z))t are martingales we 
conclude that E^M^ (,£„)] = 1 = E^M^z)]. Therefore according to Scheffe's Lemma we have 

(2.7) M t (z n ) -> M t (z) 
in L 1 (Q) as n — )• oo. 

It remains to show that f(Xf(x n )) converges to f(X^(x)) as n — > oo in Q-probability. 
Observe that together with (2.6) and (2.7) this completes the proof of 

lim P t f(x n ) = P t f(x). 

Let q t (-, •) be the density ofQ(X?(x) G • ) with respect to fi. From the strong Feller property 
of (Xt(x))t>o and another application of Scheffe's lemma one deduces that the sequence 
(qt(x n ,)) ri £N is uniformly integrable in L l (n). Hence for any e > there exists 5 > such 
that for all measurable A C S, 

fj,(A) <5=> I q t (x n ,x)n(dx) = Q(xf(x n ) G A) < e for all n and Q(xf(x) G A) < e. 
J A 

By Egorov's theorem there exists a compact set K C S such that ^{K c ) < 5 and / \ K is 
uniformly continuous, i.e. for any r/ > there is 5' > such that \z — z'\ < 5',z,z' G K => 
|/(z) - < r/- There exists iV > such that for n > N 

J^^. ®(\X?(x n ) - X?(x)\ > 5') < s 

and thus 

Q{\f(X°{x n )) - f(X?(x))\ > V )< Q{\X?(xn) - X?(x)\ > 5') + Q{X?{x n ) $ K) 

+ Q(X?(x) $K)< 3e. 

As observed earlier this completes the proof. □ 

In order to make the paper independent from the approach used by Grothaus and Klar, 
which used the machinery of generalized Dirichlet forms we add 

Proposition 2.4. If ip satisfies e~ v G L X (]R 2 ) then the measure fi introduced in (1.2) defines 
an invariant distribution for the process (X t )t>o- 
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Proof. Again we emphasize that the result is proved in [6] by different means. As mentioned 
in the introduction a direct calculation shows that the dual L of L is given by 

~l = y% ~ - + ^ • (- c :^ } ) o a + v^(o • {-_zU 

and that Lfi = 0. Setting = j^L^fi •), where by a slight abuse of notation fi is also used to 
denote the density of the measure (1.2), we get 

L» = Y% ~ cosia^ - sin(a)% + Vtffl • (~^ } ) d a . 

Completely analogous to the proof of Proposition 2.1 and Proposition 2.3 one constructs 
a conservative Feller process solving the martingale problem associated to If". From the 
conservativity we conclude by Theorem 8.5 of [20] that fi is in fact an invariant distribution. 
Observe that Theorem 8.5 in [20] - though formulated for elliptic diffusion - holds true in 
our situation with the same proof. □ 

The next result is concerned with the positivity properties of the process. 

Proposition 2.5. Let ip G C^M 2 ) and let (O, J 7 , (J r t)t>o, {Wx)xeS> {Xt)t>o) denote the unique 
weak solution of (1) on the path space Q = C([0, oo), S). If B C S is an open subset of the 
state space and if x G S then there exists T > such that for t > T 

¥ x (X t eB)>0. 

Moreover, every skeleton Markov chain is irreducible. 

Proof. The irreducibility of skeletons follows from the previous assertions together with Propo- 
sition 2.3 and section 6.1.2 in [16]. 

We first show the assertion for the case <p = 0. Let us denote by B(£,e) the ball of radius 
e > and center £. It suffices to prove that for a set B of the form B(£,e) x (q — e, a + e) 
with i G M 2 , a G §i and e > there exists T > such that for t > T 

Q x (X t eB) >0, 

where we have set x = (£, a). Assume we know that there exists a continuous path (a s ) s >o 
with values in Si such that at G (a — e, a + e) and 

(2.8) ct t G (a — e, a + e) and (^j cos(a s )ds,J sm(a s )ds^J G B(£, r) 

then there exists obviously 6 > such that for every path (a' s ) s >o with sup s \a s — a' s \ < 5 we 
also have a' t G (a — e, a + e) as well as 



J cos(a' s )ds,J sm(a' s )ds^j G B(j-,r). 



Since the support of the Wiener measure on C([0,oo),R 2 ) coincides with C([0, oo), M 2 ) we 
conclude that 

Qx(X t G B(l r)x(a-e,a + e)) > F™ & ({a' \ sup \a' s - a s \ < 5}) > 0, 

s<t 

where P^ M denotes the law of a two-dimensional Brownian motion on C([0, oo), M 2 ), starting 
at «o- Thus we only have to construct (a) s ) s >o having the properties (2.8). This can certainly 
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be deduced from control theory but we prefer a direct hands on argument. We choose the 
angle a such that the vector (cos(d), sin(d)) T is a positive multiple of the vector £ — £o ; say 
(cos(d), sin(d)) T = c(f — £o)- Let T = 1/c and T > t<i > t\ > and let (a t t 1 ' t2 )o<t<t 1 be 
defined by a/' 2 = (1 — t/t\)ao + t/t\a. For i G (£1,^2) we set c^ 1 '* 2 = a and for ti < t < T 
we set a* 1 '* 2 = (l — ^ ) a + jr* f 2 a. By choosing ii and T — £2 sufficiently small one gets 
the properties (2.8) for t = T and the extension to t > T can be done by extending a* 1 '* 2 in 
an appropriate way similiar to what has just been done. 

This establishes the assertion for cp = 0. The general case then follows from the Girsanov 
formula (2.2). □ 

We are now using ideas from the proof of Theorem 1 of [22] in order to prove the Harris 
recurrence of the fiber lay down process. Recall that (Xt)t>o is called Harris recurrent if there 
exists a finite measure 7 on the Borel cr-algebra over the state space such that for every Borel 
set A with 7(^4) > we have 



; Qf l A (X s )ds<oo^J =0 



for every x £ S. The definition of Harris recurrent for a Markov chain is analogous, the 
integral has to be replaced by a sum. {Xt)t>o is Harris recurrent if and only if the skeleton 
chain X^ = (X/ in ) n£ N is Harris recurrent for some h > 0. 

Proposition 2.6. The fiber-lay down process (Xt)t>o is Harris recurrent. 

Proof. The proof closely follows the proof of Theorem 1 in [22] and is included here for 
reasons of completeness. We first note, that the invariant distribution /x is ergodic in the 
sense, that /u is the unique invariant distribution. It follows from Proposition 2.5 that two 
different invariant distributions have to be mutually singular and Proposition 2.3 implies that 
invariant distributions have to be absolutely continuous with respect to the Lebesgue measure 
on 5. As the support of fi is the full state space the probability measure // must be unique. 
Thus the individual ergodic theorem implies that there exists a /i-null set N C S such that 
for any x £ S \ N and for any Borel set A C S 



1 " 

lim - Vp w (x,A) 
n— >oo n * — • 



n 

7=0 

Since N is also a Lebesgue null set we conclude by Proposition 2.3 that Pt(x, N) = for every 
t > and every x £ S and therefore for x £ S 

71 „ 71—1 „ 71— 1 

n- 1 Y j P h i{x,A)= / n- i y2P hi (y,A)P h (x,dy)= / n^ 1 V P hi (y, A)P h (x, dy) -> fi(A) 
7=1 J s i=0 Js\N i=0 

as n — > 00. Hence 

1 n 

lim - S2 p hi(x,A) ^ /i(A). 

77— >oo n — ' 

7=0 

for every x £ S and every Borel set A. Thus for any Borel set A with jJ,(A) > we have for 
every x £ S 

00 

^P nh (x, A) = 00, 

77 = 
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showing the recurrence of the skeleton chain. Thus according to 9.1.5 in [16] we can decompose 
the state space into disjoint sets H and N, i.e. S = H U N, such that H is a maximal Harris 
set and N is transient such that fi(N) = and H = {y S S | F v (kh < oo) = 1} and (X n h) n 
restricted to H is Harris, where kh = inf{n | X n h € H}. Since N is also a Lebesgue null 
set we conclude that P n h(x, N) = for every x £ S and n > 1. Thus, from the definition of 
H, N must be empty and (Xh n )neN Harris. Therefore the skeleton chain (Xh n )n€N is Harris 
recurrent and thus the same is true for (Xt)t>o- D 

The above results allows us to deduce the following corollary. First recall that a set C C S 
is called small for (Xt)t>o, if there exists a T > and a non-trivial measure m on the Borel 
sets of S, such that 

\fx £ C : PtO, •) > mr(0- 
The Markov transition kernel Pt is called aperiodic if for some for some small set C C S there 
exists T > such that for every x € C and t >T one has Pt(a;, C) > 

Corollary 2.7. compact subsets of the state space S are small sets for the Markov process 
(X t )t>o and the Markov process (X t )t>o is aperiodic. 

Proof. According to Propositions 2.3 and 2.6 together with Lemma 2.5 we know that (Xt)t>o 
is an irreducible strong Feller process. Thus according to Theorem 4.1 in [17] compact sets 
are petite sets and according to Proposition 4.1 in [17] petite sets are small for (Xt)t>o- The 
second assertion now follows from the fact that open sets are small and Proposition 2.5. □ 

2.2. Strong Mixing Properties. We now move on and prove ergodicity of the fiber lay 
down process, more precisely we show that it is strongly mixing in the sense that for any 
initial distribution v on S 



(2.9) lim\\F JX t G •) - ^11 = 
as well as 

(2.10) lim 

t— >oo 



%0,ao) [/(&>«*)] - / f(x)fi(dx 



S 



= 



for every / € L 2 (/i). Observe that on the one hand these assertions are stronger than the 
one in (1.3), since they do not require to look at time averages and only rely on very weak 
assumptions concerning ip. On the other hand they are much weaker than (1.3) as the do not 
provide any information concerning the rate of convergence. 

Theorem 2.8. Assume that (p satisfies the conditions of Proposition 2.3. Then the assertions 
(2.9) and (2.10) hold true. 

Proof. As is well-known strong mixing will be a rather direct consequence of Doob's theorem 
once irreducibility is established. More precisely, assertion (2.9) follows e.g. from Theorem 1 
of [22], where one has to observe that the assumption 'for a h > the transition measures 
Ph(x, •) are for x £ S equivalent' is only used in order to prove Harris recurrence for the 
skeleton chain (X n h) n£ ^ . For the fiber-lay down process (X t )t>o this has been shown in 
Proposition 2.6. Following the proof of Theorem 1 of [22] we can conclude that the skeleton 
chain (^ n / l ) n gN converges as n — > in total variation to ji and according to Theorem 3.7 in 
[21] this implies that 

lim||P x (X t G •)-m|| T v = °- 
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for ever x E S. In order to prove assertion 2.10 we rely on Korollar 3.11 in [5] (see also 
Theorem 4.5 in [3]). Since according to Propostion 2.3 we know that for every t > the 
transition kernels Pt consist of integral operators having a non-negative kernel pt(-,-). Thus 
Korollar 3.11 implies that for any / £ L p {jj,) (1 < p < oo) 



lim 

t— >oo 



Ptf 



fdfj, 



0, 



Lp(jj) 



completing the proof of the theorem. 



□ 



3. Geometric Ergodicity 

Until now it seems to be unknown, whether the rate of convergence towards the invariant 
distribution is geometric. This is at least strongly suggested by numerical experiments pre- 
sented in [6] for the cases </?(£) = |£| 2 and </?(£) = |£|. In this section we will investigate the 
question of geometric ergodicity. From the theory of Markov chains it is known that geometric 
ergodicity is strongly connected to exponential integrability of hitting times of compact sets. 
Indeed, the existence of a Lyapunov function - also called geometric drift conditions - are 
essentially equivalent to the exponential integrability of hitting times of a compact set, which 
is a small set by Corollary 2.7 (see [26] and chapter 15 in [16]). Since the angle component 
at of the process moves in a compact state space the problem consists in showing, that the 
^-component returns to compact sets fast enough. Considering the first hitting time of a ball 
in M? of radius R around it is tempting to look at the process (|&|)t>o- Assuming now that 
<p is spherically symmetric it is convenient to write the SDE (1) in polar coordinates, i.e. we 
set £ = (rcos(^),rsin(^)) and f3 = a — ip. According to formula (13) in [10] (see also p. 5 in 
[4]) this leads to the following system of stochastic differential equations 

dr t = cos(/3f) dt 

( 3 ^ dfy = (b(r t ) - —) sm((3 t ) dt + a dB t 

sin(/3t) 
dipt = , 

n 

where according to [10] b(r) = <p'(r). Equation (3.1) holds in the region, where ip is radial. 

This system has the advantage that the equation for (ipt)t>o is decoupled from the other 
( r t> /3t)t>o-P r ocess and that the equation for the radial part looks quite handy. Let us assume 
that the following assumption holds: 

A) There exists R and a constant c > such that b(r) — - > c for every r > R. 

This assumption is needed in order to make use of comparison arguments with the following 
SDE: 

dst = cos(<^) dt 

(3.2) 



d5t = csin(<^) dt + a dB t 



At this point one has to note, that we are interested in the first hitting time of the process 
(Pt, r t)t>o of compact sets of the form Si x [0,R] and therefore the singular behavior near 
will be of no concern here. 

Obviously it is not necessary for us to assume that cp is radially symmetric on the full M 2 . 
Let us thus weaken this condition slightly. We call a potential tp eventually radial, if there 
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exists a ball B(0,l) with radius I > and centre 0, such that p \ M. 2 \ B(0,l) is radially 
symmetric. 

Proposition 3.1. Assume that ip € C 2 (K 2 ) is eventually radial and assume that Assumption 
A) is satisfied. Then there exists A > such that for every (ro,/3o) € (0,oo) x [0,27r) 

i.e. the random variable tr has exponential moments of order A > 0. 

Let us for a moment assume that Proposition 3.1 holds true and let us see, how this 
implies the following theorem, which constitutes the main result of this paper. We formulate 
the theorem for the case of a gradient drift V<£> in equation (1) but we explain in Remark 3.14 
that the gradient form of the drift is not important for the following theorem. 

Theorem 3.2. Assume that ip £ C 2 (IR 2 ) is eventually radial and assume that Assumption A) 
is satisfied. Then there exists a constant v > and a strictly positive function M : S — > (0, oo) 
such that 

Vt > : ||P( &I a ) ((&>«*) e •) ~4tv ^ M{r 0: p )e- Vt 

Proof. We want to use the results from [2]. In order to do so, we use the fact that according 
to Corollary 2.7 the compact set B(0,R) x Si is small. 

Furthermore, according to Proposition 3.1 we know that a function Vr can be defined by 

" r r R( s ) 

(3.3) V R (Z,a) = l + E x=( s,a) / 

Uo 

where the 5, n > are some constants and tr(5) = inf{i > 5 \ X% G B(0,R) x §i}. Observe 
that at time 5 > the ^-component of the process travelled at most a distance 5, i.e. if £o = £ 
then | £5 — £| < 5 and therefore existence of an exponential moment for the random variable 
tr implies the existence of an exponential moment of the random variable tr{5). This shows 
that Vr is indeed well-defined. Then Vr is locally bounded and therefore by Theorem 6.2 in 
[2] the drift conditions needed in Theorem 5.2 in [2] are satisfied. Therefore we can apply 
Theorem 5.2 in [2] in order to deduce geometric ergodicity. □ 

The proof actually also allows to conclude the following slightly stronger statement. First, 
let us recall the definition of the weighted supremum spaces Bv R (S), 

B Vr (S) := { f : S ^ E | sup < 00 

I xes V R (x) 

equipped with the norm || • \\y R = sup xeS y^jj- 

Corollary 3.3. Under the conditions of Theorem 3.2 the transitions semigroup (Pt)t>o admits 
a spectral gap in the sense that for some constants C > and v > 



e xt dt 



\P t - 1 <E> /x L < Ce 



-Vt 



where 1 (8> n is the projection on the constant functions defined by 1 ® = / f dfx-X. 
Proof. This is shown in [2] to be a consequence of Theorem 3.2. □ 
In the functional analytic language of [26] and [14] this can be formulated in the following 



way 
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Corollary 3.4. Under the conditions of Theorem 3.2 the operator semigroup (Pt)t>o is quasi- 
compact in the Banach space By a (S). 

Thus it remains to prove Proposition 3.1. This will be the content of the next subsection. 
Before we start we have to emphasize once more that in this work we will be interested in the 
first hitting time of the process (Pu r t)t>o of sets of the form Si x [0, R], where R > is chosen 
large enough such that b(r) — 1/r > c for some c > and all r > R. Thus we do not need to 
worry about the 1/r-singularity of the drift in (3.1). In fact we can without mentioning this 
every time replace b(r) — 1/r by some smooth and bounded function b. 

3.1. Proof of Proposition 3.1. The proof is divided into smaller parts. 

Lemma 3.5. Let (rt,fit)t denote a solution of SDE (3.1) starting at (ro,7r) then there exists 
a reflected Brownian motion (Rt)t>o starting at such that almost surely 

|vr — fa\ < oR t 

for all t < P( 7r /2,37r/2) = inf{f > | (3 t (vr/2, 3vr/2)}. Moreover, one has almost surely 
\P pR - vr| < vr/2, where p R = mf{t > 0\R t = vr/2}. 

Proof. For the proof we set a = 1. Observe that due to 

dPt = dB t + (b(r t ) - r," 1 ) sin(&) dt 

the process (Pt)t defines a continuous semimartingale. Thus we can use the ltd - Tanaka 
formula in order to deduce that 

sgn(/3 t - vr) d/3 t + d ^^Kf^t^ 

= sgn(/3 t - vr) dB t + sgn(/3 t - vr) (b(r t ) - r^ 1 ) sin(A) dt + dL^\ 

where L^ ,7T denotes the local time of /3 in it. Observe now that for t < P( n /2,3n/2) we have 

(3.5) sgn(/3 i - vr) (b(r t ) - r^ 1 ) sin(A) < 

and that by the Levy characterization the process (J" * sgn(7r — {3 S ) dB s ^ t>Q is a Brownian 
motion. Moreover, by Theorem 22.1 in [9] we have that 

I? t ' n = - inf / sgn(/3„ - it) df3 u . 
s <* Jo 

Let us define to by 

rs rto 

(3.6) inf / sgn(/3 u - vr) dj3 u = / sgn(/3„ - vr) d(3 u 
°< s <Uo Jo 

and let t\ satisfy 

rs rti 

(3.7) inf / sgn(/3 n - vr) dB u = / sgn(/3 u - it) dB u . 
0<s<tJ Jo 

Assume that t± > to, then we by the very definition know that 

rto rto rto 



rto rto rto 

/ sgn(/3 n - vr) d/3 u = / sgn(/3 M - vr) dB u + / sgn(/3 u - ir)5 u du 
Jo Jo Jo 

rti rti rti 

< / sgn(/3 n - vr) d/3 u = / sgn(/3 u - vr) dB u + I sgn(/3 M - tt)S u du, 
Jo Jo Jo 



(3.8) 

where 5 U := (b(r t ) - r 4 _1 ) sm(/3 t ). 
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We might assume to have a strict inequality, since otherwise we could take to = t±. Now a 
strict inequality in (3.8) implies that 

ft\ rt± 

/ sgn(/3„ - tt) dB u >- sgn(/3 u - tt)5 u dy > 0. 

J to J to 

But this gives us 

"tl rk) rti 



rti fto rti 

/ sgn(/3 n - 7r) dB u = / sgn(/3„ - tt) dB u + / sgn(/3 M - tt) dB u 
Jo Jo Jt 

/•to 

> / sgn(f3 u - tt) dB u 
Jo 



contradicting our choice of t\ in (3.7) and the assumption t\ > to. 
In the case t\ < to we proceed as follows : We have 



ft ft fS 

/ sgn(/3 M - it) dB u + / sgn(/3 u - tt)6 u du - inf sgn(/3 u - tt) df3 u 

Jo Jo s ^Jo 

rt ft fto 

I sgn(/3 n - 7r) dB u + / sgn(/3 n - 7r)(5 u du - I sgn(/3 u - 7r) 

io io Jo 



fto 

+ / sgn(/3 u - ir)5 u du 
Jo 

ft ft fS 

< / sgn(/3 u - tt) dB u + / sgn(fi u - 7r)5 u du - inf sgn(/3 u - 7r) dB u 
Jo J to s ^Jo 

ft fto 

< / sgn(/3 n - tt) dB u - I sgn(/3 u - 7r) 
io io 

where we used (3.5) in the last step. This gives in the case t± < to 

(3.9) \fo - tt\ < [ sgn(/3 u - tt) dB u - inf / sgn(/3 u - tt) dB u , 

Jo s <*io 

and since according to Theorem 2.34 in [13] the right hand side is distributed as the reflected 

Brownian motion we arrive at the assertion of the lemma. □ 

In a completely analogous way one proves 

Lemma 3.6. Let (rt,Pt)t denote a solution of SDE (3.1) starting at (ro,2-7r) then there exists 
a reflected Brownian motion (Rt)t>o starting at such that almost surely 

\2ir > oRt 

for all t < P( 7r /2,3 7r /2) = inf{t > | fit £ (tt/2,3tt/2)} . Moreover, one has almost surely 
I/Vap W2 ,3./ 2) -*\> V2, where p R = m£{t > 0\R t = tt/2}. 

3.2. Special Case. Assume that b(r) = c+l/r for r > R. Then the process (Pt)t>o satisfies 
dfit = o~dB t + csin(/3 t ) dt up to the random time 

J t r = inf jt > | J cos(/3 t ) dt< r\. 

This is a very special case, but it already contains all features necessary to understand the 
general situation. Moreover, we make use of the process (7t)t>0) which is defined in the 
following way. 
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i) 7o = tt; It behaves like a Brownian motion reflected at tt up to the first hitting time 
of 3vr/2. 

ii) Then jt behaves like the solution of dfy = crdBt + c sin(/3t) dt until it hits tt or 2tt. 
hi) If in step ii) 7$ hit tt, then we proceed as in item i). If in step ii) jt hit 2tt, then we the 

process jt is continued as a Brownian motion reflected at 2tt up to the first hitting 
time of 37r/2. After this time we proceed as in item ii). 

Let us denote by (f3t)t>o the solution of dj3 t = adB t + c sm(/3 t )dt started at tt and we set 
ft = J*g cos (P s ) ds. The process 0t)t>o is still considered as a process on the circle, i.e. we 
identify and 2ir. Thus up to the first hitting time tr = tr of the set §>i x [0, R] the two 
processes (f3 t , r t )t>o and (fit, ^)t>o coincide if /?o = /3o and tq = tq. We now define recursively 
the stopping times 

a = = a , o x = mf{t > | ft £ (0, 2vr)}, a 1 = mf{t > | ft = vr}, a 1 = a l + a 1 o 6 ai 
and 

The analog stopping times are defined for the process (7)t>o and by a small abuse of notation 
we denote these also by (a l )i e ^ and ((7i)i^ . 

We say that the process (/3t)t>o (resp. (7t)i>o) started at tt completes a cycle if it returns 
to tt after having visited 2tt. Thus the stopping time a 1 describe the time of completion of 
the i-th cycle. 

Moreover, for some interval we denote by Aq the smallest Dirichlet-eigenvalue of the oper- 
ator — \-^2 — csin(x)^, which can also be characterized by 

Aq = - lim \ log sup ¥ z (pi > t) , 

t^oo t ze i 

where pj = inf{t > | p\ ^ /} denotes the first exit time. We set 

\ _ rn ; T1 /'\(0,T) x (7r/2,37r/2) y(*,2ir)\ 

A — mm(A , A , A J. 

Lemma 3.7. The random variables 



Xi = cos(P s ) ds, i = 1,2, .. . 



form an i.i.d sequence of random variables having the following properties: 

• The random variables Xi (i = 1,2, . . . ) have a negative expectation, i.e. 

• For all A < Ao we have 

% 0= Je^l] < oo. 

The same properties hold true for the sequence (X?)^ where X] = X^-i cos(7 s ) ds, i = 
1,2,... 

Remark 3.8. Though we are not going to present a proof here it seems possible to show, that 

Es [AM 
lim ^° 1J = 0. 

o-^oo His _ \o~ 
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Heuristically this is to be expected since for a large diffusion constant the influence of the drift 
becomes eventually negligible. Such an assertion might be finally used in order to argue that 
enlarging the diffusion constant a does not lead to faster rate of convergence. Instead one 
even has to expect that the rate of convergence will become slower. This is due to the fact 
that the value E^ Q=7r [Xj] controls in a certain sense the 'speed' with which the process (Xt)t>o 

returns to a ball around the origin. Bigger values o/ E^ Q=7r heuristically lead to longer 
return times to the ball and thus one expects to get a smaller spectral gap. 

Proof. The fact that the sequence is i.i.d is a direct consequence of the strong Markov 

property of the process (f3t)t>o- 

The first assertion can be shown in the following way. We first observe that 



(3.10) 



cos(/3 s ) ds, 



E 



/9o=t 



cos(/3 s ) ds, 



+ E 



Po=2tt 



cos(/3 s ) ds, 



We will now compare the first and the second summand in equation (3.10). By the strong 
Markov property we have 



E 



A)=2tt 



cos (/?<;) ds 



and then 



P0=2n 



E 



A)=27T 



P(37r/2,57r/2) 



cos(/3 s ) ds 



o 



+ E 



A) = 27T 



cos(/3 s ) ds 



E 



/3 =2tt 



P(37r/2,5ir/2) 



cos(/3 s ) ds 



o 



+ 2 %o=3tt/2 



cos((3 s ) ds 



2 po=t/2 



P(37r/2,57r/2) 



cos(/?s) ds 



cos((3 s ) ds 



E 



/3o=2tt 



I 



P(3ir/2,57r/2) 



cos(/3 s ) cis 



+ 2 E /3o=3tt/2 



P(7T,27r) 



COS 



(A) da] 



(3.11) 



+ ^ P /3o=3/2.(^ 1 > P(7r,2 7 r))E 27 r 



cos (/3s) 



+ -Es /0 
2 po=t/2 



P(0,7r) 



cos(/3 s ) <is 



1 



+ o P /3o=V2(^ >P(0,.))E 27r 



cos(/3 s ) ds 



where pi denotes the first exit time from the interval I. Thus with d\ = 1 — Eg o=37] y 2 (cx 1 > 



P(o,tt)) we get using obvious symmetry properties 



(3.12) 



E A,=2. 



cos(/3 s ) (is 



J /3 =27T 



P(3ir/2,5ir/2) 



cos(/3 s ) (is 



o 



+ E 



/3o=37r/2 



P(7T,27r) 



cos(/3 s ) ds 
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Analogously we get 
(3.13) 

cos(/S s ) ds 

FP(n /2,3tt/2) 



Es 

P0=1~ 



Es 

PO = TT 



E 



/3 =7T 



2 po=7i"/2 



Es 

P0=T 



cos (/3s) ds 



P(7r/2,37r/2) 



cos(/3 s ) ds 



+ E 



/9o=t 



cos (/?<;) (is 



P(ir/2,37r/2) 



cos(/3 s ) (is 

P(ir/2,3?r/2) 



+ 



cos(/3 s ) (is 



Uo 



cos((3 s ) ds 







+ -Es /0 

2 P0=37r/2 



P(0,7T) 



cos(/3 s ) (is 

P(7T,2ir) 

cos(/3<j) ds 



1 



+ /^^(^ > P(0,,r))E* 



CI 



cos(/3 s ) (is 



1 



+ 2^0=3^/2(^1 > P {n ,2n))^0 Q=v 



cos(/3 s ) ds 



and with d 2 = 1 - ^ = 37r / 2 > P(ir^r)) 



E 



(3.14) 



/3o=7r 



cos(/3 s ) ds, 



do 1 ! E 



2 l^/3 =vr 



E /3 =3tt/2 



P(7r/2,37r/2) 



cos(/3 s ) ds 





P(7T,27T) 



cos(/3 s ) ds 



Using Lemma 3.5 together with the properties of the cosine we conclude that 



(3.15) 



E 



/3 =7T 



P(7r/2,3ir/2) 



cos(/3 s ) ds 



> -Es 

— P0=7T 



P(7r/2,37r/2) 



cos(-B s ) ds 



where (£><,) denotes a Brownian motion in [0, 2ir] with periodic boundary conditions and that 



E 



I3 =2n 



P(37r/2,57r/2) 



cos(/3 s ) ds 



< E A,=2. 



P(ir/2,37r/2) 



cos(£? s ) ds 



Therefore 



E 



ft> = 7T 



P(vr/2,37r/2) 



cos(/3 s ) ds 



< E 



/3o=2tt 



P(3ir/2,57r/2) 



cos(/3 s ) ds 



Using Theorem 3.10 in the Appendix or alternatively explicit expression for hitting probabil- 
ities from Theorem 23.7 in [9] it is straightforward to deduce that 

P /3 =37r/2( a ' i > P(*,2*)) < IP / 3 =37r/2( (T l > P(^)) 

and therefore d^ 1 < d^ 1 . 

Now comparison of (3.12) and (3.14) using (3.15) as well as the following Lemma 3.9 to- 
gether with d^ 1 < dg -1 allows to deduce the first assertion. 

In order to prove the second assertion observe that |ATi| < a 1 . Thus we get using the strong 
Markov property again 



(3.16) 



According to our choice of Aq the right hand side remains finite for A < Aq. 
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The assertions for (X^)^ =1 can be proved in an analogous way. The first assertion remains 
true without any change. For the second assertion observe that has exponential moments 
up to Ao since for every / = (0, tt), (tt/2, 3ir/2), (tt, 2ir) it follows in a straightforward way 
from Theorem 3.15 -recalled in the Appendix - that 



- lim - log sup ({7, E I V0 < s < t}) > A 

t^oo t x&1 



giving the required estimates. 

Lemma 3.9. We have 
a) 

b) 



□ 



E 



3n/2 



P(vr,2ir) 



cos(/3 s ) ds 



E 



tt/2 



P(0,7r) 



cos(/3 s ) ds 



< 0. 



< 0. 



Both assertions remain obviously true, if (fit)t>o is replaced by ("ft)t>o- 

Proof. We only prove assertion a) and in order to simplify the formulas we set a = 1 and 
c = 1. In order to prove the general case only notational changes are necessary. We have 



E 



(3.17) 



3tt/2 



cos(/3 s ) ds 



E 



3tt/2 
oo 



cos(/? s )l {p( ^ )>s} 



ds 



o 



E 37r/2 [cos(&)l 



{P(tt,2tv)>S 



}] dS - 



Thus it is sufficient to show that E 37r / 2 [ cos (/5t)l{p (?r 97r) ><}] < for t > 0. By Girsanov's 

theorem we have 

(3.18) 

13V2[COS(A)1 W) > ] = E3 V2 [cOs(^) e /o^)^-U t sin 2 (B fl )„ 1{ ^^ >t}] 
= ^08(377/2)^^ ^ co ^ Bt ^cos(B t ) e -y^co S (B 3 )ds~lJ [ i ; S mHB s )ds 1{ ^ 

where we used, that according to the Ito formula 

- cos(B t ) = - cos(B ) + I sin(i? s ) dB s + - J cos(B s ) ds. 

Jo 2 Jo 

Finally by considering the excursion straddling t and observing that the cosine is negative in 
[tt,3tt/2) and positive in [3tt/2,2it) we conclude that 

(3.19) 

E3 w/2 [cOs(S t ) e - COS ^) e -|/o t cos(B s )d S -i / ^i^(B a )d, 1{ ^^ >t}] 

= E 3n/2 [cos(B t )e- co ^e-y» cos ^ ds -yo si » 2 ^ ds l {p( ^ >t} ;B t < 3tt/2] 
J +E 3n/2 [co S (B t )e- cos ^e-lfo^mds-y^in^ > ^/ 2 ] 

= h + I 2 < 0. 

In order to see this note first that we are using the symmetry of the expression sin 2 (i? s ) ds 
in order to restrict ourselves to the expressions containing the cosine. Then one should 
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observe that the contribution of the path (B s ) s< ^ t -j with l(t) = sup{u < t \ B u = 3n/2} to 
the integrands in I\ and I2 in (3.19) is the same and that therefore only the contribution of 
{B s )u t \ <s<t makes the difference. This completes the proof of the assertions. □ 

Define (S n ) n >i to be the random walk corresponding to the i.i.d sequence i.e. S n = 

Y^i=\Xi. Moreover let T_ denote the time 

T_ = inf{n > I S n < R}. 

Similiarly, we define (S n ) n >i to be the random walk defined by Sn = Y17=i ^7 an< ^ ^ ^1 
be the first hitting time of [— 00, R]. We first note that our results imply the existence of 
exponential moments for T_ and TZ, respectively. 

Lemma 3.10. Let v denote an initial distribution on [R,oo] having exponentially decaying 
tails, i.e. there exists cq > and c\ £ (0,1) such that u((t,oo)) < cqc\. Then there exists a 
strictly positive X' depending only on c± such that 

^u[e X ' T -] < 00 and E^[e A ' T -] < 00 

Proof. By shifting the problem using translation invariance of the random walk we can assume 
that R = 0. First we note that according to [15] we have for every initial point sq > R = 

E S0 [e XT -} <oo 

for some A > 0. Let us fix such an initial point sq and A > A' > 0. Since the measure v is 
supposed to have exponentially decaying tails we have that 

V/c > : u((ks , (k + l)s )) < cqc\. 

Define the times Tk SQ to be 

T kso = infjra > | S n < ks } 

and observe that when started at (k + l)so the hitting time of the negative real numbers can 
be bounded by first waiting for Tfc so , then starting at ksQ we wait for ?(fc_i) So and so on. This 
gives 

Therefore 

n 
k=0 

if A' > is taken to be small enough. Observe that the choice of A' depends only on c\. The 
same arguments obviously show the assertion concerning T]_. □ 

We denote the stochastic process ( / O *cos(/3 S ) ds ) by (f t ) and we set 

r R = mf{t>0\r t <R}. 
Moreover, r[ denotes Jq cos(7 s ) ds and is defined to be the first hitting time 

tI = mi{t > I rj < R). 
The assertion of Proposition 3.1 for the simple case thus follows from the next Lemma. 
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Lemma 3.11. Let v denote some initial distribution on [R, oo) satisfying v((t, oo)) < cqc\ 
for some cq > and c\ S (0, 1). Then for some A > depending only on c\ and we have 

E„ [e Affl ] < oo and E u [e Xr R] < oo 

Proof. According to the construction of the random walk (5 n ) n£ N we conclude that 

T\-aa,R] = n implies tr happened before completing cycle n. 

whenever the random walk S n and the process rt are started from the same initial distribution. 
This means that 

oo 

(3.20) E u [e XrR ] = ^P„(T_ = n)E u [e XTR ] n . 

n=l 

Using the fact that according to Lemma 3.10 the sequence (F U (T_ = n)) ne pj decays exponen- 
tially we can choose A > such that the right hand side of (3.20) remains finite. The same 
proof applies to the second assertion. □ 

Corollary 3.12. Assume that outside a compact set ip is radial and that outside this compact 
set b{r) = c + 1/r. Then the assertion of Theorem 3.2 holds true. 

Proof. It remains to show that suitable sufficiently large R > and every (ro, /3q) € (0, oo) x Si 
we have E( ro [e TR ] < oo. Observe that by the strong Markov property we have 

E (ro,/3 ) [e AT *] <E K ro)W) [e A ^], 

where v/3 0l r is some distribution on (it!, oo). Moreover, as the tails of vp Q ,r can be dominated 
by the tails of the distribution of the first hitting time of ir by the process (A)t>o we conclude 
that there exists a constant Co = co i( 3 , ro > (possibly depending on (3q and ro) and a constant 
ci £ (0, 1), which does not depend on (3q and ro such that for I > R 

Therefore there exists some A > such that for every pair (/3o 5 r o) we have E( rQ)/ g ) [e Arfl ] < oo, 
showing that Proposition 3.1 is true in this special situation and thus Theorem 3.2 as well. □ 

Remark 3.13. The analysis of this special case demonstrates the dynamical properties of the 
fiber lay down process in a rather clean way. In order to extract the net- 'drift' towards of 
the radial component one has to wait until a full cycle is completed. Then each completed 
cycle can be used as an i.i.d. increment in a random walk which has then negative drift. The 
identification of the possible translation of the diffusion hitting time problem to a random walk 
problem constitutes the essential insight. 

3.3. General case. Recall that by definition the process {^ftjt starts from ir as a reflected 
Brownian motion until it reaches ir/2 or 37r/2, then it behaves as the solution of dXt = 
dP>t + csin(X s ) ds =: dBt + 8{X s )ds until it hits ir or = 2ir, from there it starts as a 
reflected Brownian motion at 2tt until it hits tt and repeats this behavior. We have seen, that 
with rj = J * cos(7 s ) ds and 

tI = inf{t > | rj < R} 

we have 

- lim -logPJrZ > t) > 

<->oo t ' 
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for every initial distribution v on (0, oo) having exponential moments. The main goal in this 
subsection consists in proving that it is possible to compare the general case with the process 



3.3.1. A Coupling Construction. We are now going to construct a process (jt)t>o i n a way, 
which allows a comparison to the process (f3t)t>o- More precisely, ("ft)t is constructed in such 
a way that 



We want to construct the coupling for one complete cycle of (jt)t, i-e. we want to construct 
such a coupling up to the time a 1 . The idea is easy to grasp, but as often with coupling 
constructions, the formal presentation tends to hide the simple underlying intuition. The 
idea consists essentially in a comparison of the process (5 with a process which has a weaker 
drift to tt and which thus will have a bigger distance to tt. Then one can conclude using 
the properties of the cosine that the radial process r will be smaller than the corresponding 
'radial component' of the comparison process. 

The idea is the following: We use one of the following rules. 

• If the processes ("yt)t>o and (fit)t>o are at a given time both at the point tt or at 2tt 
then we use Lemma 3.5 and Lemma 3.6, respectively in order to make sure that the 
distance of the 7-process to tt dominated the distance of the /3-process to tt. 

• If the /3-process is at the point tt and the 7-process is in (tt, 2tt] we wait until the 
distance of the 7-process to tt coincides with the distance of the /3-process to tt. At 
this time either both processes are in tt or one is in the situation of the following 
item. The same procedure is applied in the case, where the 7-process is at 2tt and the 
distance of the /3-process to tt is strictly smaller than tt. 

• If at a time r the 7-process is in the interior of (tt, 2tt) and the /3-process is in (0, tt) U 
(tt, 2tt) with |/3 T — tt\ < \ j T — tt\ then one can use Theorem 3.15 in order to conclude 
that distance of the 7-process to tt remains dominated the distance of the /3-process 
to tt, at least until one hits the boundary and repeats the procedure outlined in the 
second item. If the processes 7 and /3 are both in (tt, 2tt) then the application of 
Theorem 3.15 is direct, if the process /3 is in (0,tt) and the process 7 is in (tt,2tt) 
then one can use the symmetry properties of the drift of the 7-process in order to 
conclude necessary domination of distance of the /3-process to tt by the distance of 
the 7-process to tt via Theorem 3.15. 

Let us give now the first steps in the construction on a more formal level. Let v be an initial 
distribution on [R, 00) having exponentially decaying tails and let £ ~ v be a random variable 
which is distributed according to v. 

a) We start our process (flt,r t )t with initial conditions r$ = £ and (3q = tt. According to 
Lemma 3.5 there is a Brownian motion (Rt)t>o started in tt and reflected at tt such 
that the distance of /3f to it is smaller than Rt and we wait up to time p\ when the 
process (7t)j>o = (Rt)t>o hits the point 3tt/2. 

b) For t > pi we proceed as follows. On the probability space supporting the process 
( r t,flt)t>o we can solve the stochastic differential equation drft = dB t + cs\ii(^ t ) dt, 
where (Bt)t>o is the Brownian motion driving the process (rt,fit)t>o- Observe that 
the solution (7t)t>o is pathwise unique. We wait up to the random time P2 = P2,i^P2,2, 
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where 

p 2 ,i = inf{i | 7t = 27t} and p 2 ,2 = inf{i | 7* - 7r = |/3 t - 7r|}. 

Observe that up to the stopping time P2 the distance of /? to tt does not exceed the 
distance of 7 to n. 

c) If />2 — p2,i then we extend (74)4 as a Brownian motion )t>o at 27r which is reflected 
at 2tt and wait until p 3j i = pz,\,i A P3,i,2 ; where 

= inf{i I = 7r} and p 3t i >2 = inf{t | i?J - tt = \/3 t - 7r|}. 

If P2 = P2,2 then we extend the comparison process (74 )t as a reflecting Brownian 
motion i?j if 7 P2 2 _ £ {0, 27r} or let it behave as a solution to djt = dB t + csin(7 t ) dt, 
and we wait up to the time p^^ = ^3,2,1 A /5 3 ,2,2, where 

P3,2,l = inf{t I 7t G {0, 27r} or R 2 t = tt} and 

P3,2,2 = hif{t I -Rj - TT = |/3 t - 7r| (7 t - 7T = |/3 t - 7r|)} 

meaning that we wait if we hit 7r in the case that we move according to a reflected 
Brownian motion or 0, 2ir otherwise or alternatively until the comparison process has 
the same distance to it as the process /3 . 

Observe once again that due to Theorem 3.15 we can conclude that the distance of 
the /3-process to it does not exceed the distance of the 7-process to it. 

Let us set p 3 = ^3,1 A p^ 2 - 

d) At time p 3 there are three possibilities: If the /3-process is at the boundary and the 
7-process is in the interior, then we wait until either the distance of both processes to 
tt coincides or if the 7 hits the boundary. If both processes are at the boundary we 
can make use of our Lemmas and wait until the 7-process hits tt and if both have the 
same distance to tt we can use Theorem 3.15 until one process hits the boundary. 

In all cases the distance of the 7-process to tt dominates the distance of the /3-process 
to tt. Then this procedure is repeated in the same way. 

We emphasize, that the items a) - d) do not allow to construct the coupling for all times. 
Thus we have to observe that by construction of the coupling the process (7t)t>o is always 
between the process ((3t)t>o and the point 2tt and we know, that the time a 1 is finite for 7^ 
and 7 CT i = tt. Therefore at time a 1 we have j3 a i = 7^1 = tt and we can start again from item 
a) replacing tq by r a i . 

Iteratively we obtain two processes: the process (r t ,f3 t )t>o and (7t)t>o such that for all 
t > we have 

Vt > : \tt- f3 t \ < k-7i|. 

Using this together with the fact that [tt, 2tt] 3 I 1— > cos(| Z| ) is increasing it then becomes clear 
that 

rj > r t . 

Finally we conclude that 

t r = inf{t >0\r t <R}<Tl = inf{t > | r t 7 < R} 

and since according to Lemma 3.11 the random variable has exponentially decaying tails 
the proof Proposition 3.1 is completed. 
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Remark 3.14. We would like to end with the following observation : The gradient form 
of the drift in (1) was dictated by the aim to have an explicit candidate for the invariant 
distribution. The apriori knowledge of the existence of an invariant distribution was used 
in Theorem 2.8 as an essential ingredient. We want to emphasize that this is not true for 
Theorem 3.2. For fixed T > the analog of the function Vr, defined in (3.3), satisfies (see 
Theorem 6.2 in [2]) 

P s Vr < \{s)V R + bl- m ^ )x[QM : 

where X(s) is bounded for s G (0,T] with X(T) < 1 and b £ (0,oo). Therefore by Theorem 2.1 
of [2] we can conclude the existence of an invariant distribution for the fiber lay down process. 

Thus one might hope that the approach presented in this work is also applicable to the case 
of a moving conveyer belt. A detailed investigation of this is deferred to the future. 

Appendix 

In this section we recall Theorem 1.1 from chapter VI [8]. Suppose we have two real 
continuous functions b± and b 2 defined on [0, 00) x K such that 

(3.21) Vt > 0, x <E M : h(t, x) < b 2 {t, x) 

and let (fi, J-, P) be a probability measure with a reference familiy (J-t)t>o- 

Theorem 3.15. Suppose that we are given the following processes: 

1) two real (T)t>o- adapted continuous processes x\{t,u) and x 2 {t,u) 

2) a one- dimensional {Ft)t>o-Brownian motion B{t,uj) such that Bq = a.s.. 

3) two real (J r t)t>o- adapted well-measurable processes (3\{t,uj) and ^{t^uj). 
We assume that they satisfy the following conditions with probability one: 

(3.22) Xi(t)-Xi(0) = B(t)+ [ Pi(a)ds,i = 1,2, 

J 

(3.23) xi (0) < s 2 (0) 

(3.24) Pi(t) < &i(t,xi(*)) for every t > 

(3.25) p 2 (t) > b 2 {t, x 2 {t)) for every t > 0. 
Then with probability one, we have 

(3.26) X! (t) < x 2 {t) for every t > 0. 

Furthermore, if the pathwise uniqueness holds for at least one of the following stochastic 
differential equations 

(3.27) dX{t) = dB{t) + bi{t, X(t)) dt, i = 1,2, 

then the same conclusion (3.26) holds if (3.21) is replaced by the weakened condition 

(3.28) h(t, x) < b 2 (t, x) fort>0,x e R. 
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